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Abstract 

In these lectures, a variety of non-equilibrium transport phenom- 
ena are introduced that all involve, in some way, elastic manifolds be- 
ing driven through random media. A simple class of models is studied 
focussing on the behavior near to the critical "depinning" force above 
which persistent motion occurs in these systems. A simple mean field 
theory and a "toy" model of "avalanche" processes are analyzed and 
used to motivate the general scaling picture found in recent renormal- 
ization group studies. The general ideas and results are then applied 
to various systems: sliding charge density waves, critical current be- 
havior of vortices in superconductors, dynamics of cracks, and simple 
models of a geological fault. The roles of thermal fluctuations, defects, 
inertia, and elastic wave propagation are all discussed briefly. 



I. Introduction 

Many phenomena in nature involve transport of material or some other quan- 
tity from one region of space to another. In some cases transport occurs in 
systems that are close to equilibrium with the transport representing only 
a small perturbation such as flow or electrical current in a metal, while in 
other cases it involves systems that are far from equilibrium such as a land- 
slide down a mountain, or a drop of water sliding down an irregular surface. 
Sometimes, particles or other constituents move relatively independently of 
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each other like the electrons in a metal, while in other situations the inter- 
actions play an important role, as in the landslide and the water droplet. 

If the interactions are strong enough, all the particles (or other con- 
stituents) move together and the macroscopic dynamics involves only a small 
number of degrees of freedom. This is the case for a small water drop on, e.g. 
wax paper, which slides around while retaining its shape. But if the interac- 
tions are not so strong relative to the other forces acting on the constituents, 
then the transport involves in an essential way many interacting degrees of 
freedom. This is the case for a larger water drop on an irregular surface for 
which the contact line between the droplet and the surface continually de- 
forms and adjusts its shape in response to the competition between the surface 
tension of the water and the interactions with the substrate QT5|| . Such 
a moving drop and the landslide are examples of non-equilibrium collective 
transport phenomena, which will be the general subject of these lectures. 

This is, of course, an impossibly broad subject! We must thus narrow 
the scope drastically. Although the range of systems discussed here will, 
nevertheless, be reasonably broad, we will primarily focus on systems in 
which the interactions are strong enough so that the transported object (or 
at least some part of it) is elastic. We will use this in a general and somewhat 
loose sense that the transported object has enough integrity that if one part 
of it moves a long distance then so, eventually, must the other parts as well. 
Thus the fluid drop is elastic if it does not break up — i.e. its perimeter retains 
its integrity — while a landslide is not elastic as some rocks will fall much 
further than others and the relative positions of the rocks will be completely 
jumbled by the landslide. 

We will be interested in systems in which the medium in which the 
transport occurs has static random heterogeneities ( "quenched randomness" ) 
which exert forces on the transported object that depend on where it is in 
space. 

Examples we will discuss are: interfaces between two phases in random 
media p8|P], such as between two fluids in a porous medium or do- 



main walls in a random ferromagnetic alloy; lattices of vortices in dirty type 
II superconductors |34]]; charge density waves which are spatially periodic 
modulations of the electron density that occur in certain solids {| Q ; and 
the motion of geological faults In addition to the contact line of 



the fluid drop already mentioned |?6],|16j, another well known — but poorly 



understood — example that we will, however, not discuss is solid-on-solid fric- 
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tion. 

In all of these systems, a driving force, call it F, can be applied which 
acts to try to make the object move, but this will be resisted by the random 
"pinning" forces exerted by the medium or substrate. The primary questions 
of interest will involve the response of the system to such an applied driving 



force |47|| . If F is small, then one might guess that it will not be sufficient 
to overcome the resistance of the pinning forces; sections of the object would 
just move a bit and it would deform in response to F, but would afterwords 
be at rest. [Note that in most of what follows, we will ignore fluctuations so 
that the motion is deterministic and the objects can be said to be stationary] 
If the force is increased, some segment might go unstable and move only to 
be stopped by stronger pinning regions or neighboring segments. But for 
large enough F, it should be possible to overcome the pinning forces — unless 
they are so strong that the object is broken up, an issue we will return to 
at the end — and the object will move, perhaps attaining some steady state 
velocity v. 

Basic questions one might ask are: is there a unique, history independent 
force, F c separating the static from the moving regimes? How does v depend 
on F (and possibly on history)? Are there some kinds of non-equilibrium 
critical phenomena when v is small? How does the system respond to an 
additional time or space dependent applied force? These are all macroscopic 
properties of the system. 

But we will also be interested in some microscopic properties: how can one 
characterize (statistically) the deformations of the object when it is station- 
ary ||38|| ? The dynamic deformations and local velocities when it is moving? 



The response to a small local perturbation? etc. 

Motivated by possible analogies with equilibrium phase transitions [|33 
we can ask if there are scaling laws that might obtain near a critical force 
which relate, for example, the characteristic length scale L, for some process, 
to its characteristic time scale, r, via a power law relation of the form: 

t ~ U (1) 

Trying to answer some of these questions — and to pose other more pointed 
questions — is the main aim of these lectures. In the next few sections a 
particular system and its natural (theoretical) generalizations will be studied 
and tools and ideas developed. In the last section, these are tentatively 
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applied to various physical systems and some of the complicating features 
left out of the simple model systems are discussed. This leads naturally to 
many open questions. 



II. Interfaces and Models 

In order to develop some of the general ideas — both conceptual and compu- 
tational — we will focus initially on an interface between two phases that is 
driven by an applied force through an inhomogenous medium p8|,P|. The 
essential ingredients of a model of this system are: the forces of sections of 
the interface on nearby sections, i.e. the elasticity of the interface caused by 
its interfacial tension; the preference of the interface for some regions of the 
system over others due to the random heterogeneities; and some dynamical 
law which governs the time evolution of the local interface position. 

We will initially make several simplifying approximations, which we will 
come back and examine later. First, we assume that the interface is not too 
distorted away from a flat surface normal to the direction (z) of the driving 
force so that its configuration can be represented by its displacement field 
u(f) away from a flat reference surface. The coordinates R = (x, y, z) of 
points on the surface are then 

(x, y) — f and z — u(r). (2) 

Second, we will assume that the dynamics are purely dissipative i.e. that iner- 
tia is negligible — a good approximation in many physical situations. Keeping 
only the lowest order terms in deviations from flat, we then have 

V ^&T = F + t] ~ u [f > u{f > t)] (3) 

with F representing the driving force on the interface, f p (R) representing 
the random "pinning" forces of the heterogeneous medium on the interface 
which we assume for the present are not history or velocity dependent; rj a 
dissipative coefficient; and 

a(r, t) = J dr J* dt'J(r -r',t- t') [u{f', t') - u(r, £)] (4) 
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interface 



Figure 1: Schematic of a one-dimensional interface in a two dimensional 
disordered system illustrating the forces acting on the interface. 



the "stress" on the interface from its elasticity which is "transmitted" by the 
kernel J(f,t). Short range elasticity of the interface corresponds to 

Ja5{t)V 2 5{f). (5) 

A schematic of such an interface and the forces acting on it is shown in Fig |T[ 
Keeping in mind some of the other problems of interest f36f,f6|, |30| in 



addition to interfaces, we will abstract to a more general problem of a d- 
dimensional elastic "manifold" — d = 2 for the interface — with more general 
interactions, which can be long-range, embodied in J(r,t). In addition to 
the form of J(r,t), the system will be characterized by the statistics of the 
pinning forces which impede interface motion near points where the interface 
has lower (free) energy; f p (R) will generally have only short-range correla- 
tions in space, i.e., in both u — u' and r — f. Even with these simplifying 
assumptions, the model Eqs £|) is impossible to analyze fully due to the 
non-linearities implicit in the u dependence of f p (f,u). Nevertheless a lot of 
the qualitative behavior can be guessed. 

If the driving force is sufficiently small, then it will be insufficient to 
overcome the pinning forces. But if F is increased slowly, it may overcome 
the pinning of some small segment of the interface which can then jump 
forwards only to be stopped by stronger pinning forces or by the elastic 
forces from neighboring still-pinned parts of the interface. But if the drive 
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is larger, the neighboring regions may themselves not be strongly enough 
pinned to resist the increase in stress from the jumping section and may 
themselves jump forward leading to an "avalanche" of some larger region of 



the interface; this process might or might not eventually stop ||38| . If the 
force is large enough — and certainly if it exceeds the maximum f p — then it is 
not possible for the interface to be pinned and the interface will move forward 
with some average velocity v, albeit via very jerky motion both in space and 
time. 

In addition to the basic questions raised in the Introduction, we will 
be interested, in the regime where avalanches can stop, in the statistical 
properties of the sizes and dynamics of the avalanches [[38] || || . 

The existence of a unique critical force can be established by a simple 
convexity argument that is valid if J(r, t) is non- negative . If two con- 



figurations u a and Uj, of the interface at time to have the property that one 
is "ahead" of the other, i.e. u a (f,t ) > Ub(r,t ) for all r; then u a will be 
ahead of Ub at all later times. This can be seen by assuming the contrary 
and then considering the putative first time t\ > to at which there is a point 
of contact; say u a (fi,tx) = «&(fi,ti) at some r\. Then 

d 

— [u a (fi,t) - u b {fi,t)] \ t=tl = a[n, {u a }] - a[n, {u b }] 
= J dftpdt' J(r 1 -f*,t 1 -t , )[u a (r , ,t , )-u b (r , ,t')} > (6) 

since the pinning force at f*i is the same in both configurations and therefore 
cancels out. By assumption the last expression in Eq is positive as long 
as J is non- negative so that for t > ti, u a is again ahead of Ub violating the 
assumption. 

The condition that 

J(f,t)>0 (7) 

for all r, t plays an important role in the theoretical analysis and frequently 
also in the physics of these types of systems. We will refer to models with 
this convexity property as monotonia, they have the property that if the 
displacements and the driving force F(t) increase monotonically with time, 
then so will the total "pulling force" — see later — on any segment. Except in 
the final section we will focus solely on monotonic models. 

We have shown that in monotonic models one configuration that is ini- 
tially behind cannot "pass" another that is ahead of it p6f; therefore station- 
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ary and continually moving solutions cannot coexist at the same F; therefore 
F c is unique. This is a big simplification and one that will not occur gen- 
erally, in particular not in some of the systems that we discuss in the last 
section. 

For forces well above F c , one can use perturbative methods to study the 
effects of the random pinning and compute, for example: the mean velocity, 
v(F) , the spatio-temporal correlations of the local velocities, and responses to 
additional applied forces . Near F c , however, life is much more complicated 
as is usually the case near critical points — but for conventional equilibrium 
critical points the theoretical framework for dealing with the complexities 
is well established ||33|| . The most interesting behavior occurs in the critical 



regime; in particular one might expect processes — such as avalanches — to 
occur on a wide range of length and time scales. In order to make real 
progress, we will have to — at least initially — make further simplifications or 
approximations. 

One of the lessons from equilibrium critical phenomena is that analyzing 
simple models exactly or by controlled approximations is more useful than 
analyzing more realistic models by uncontrolled approximations; we will thus 
take the former route. But we must first find some clues as to what simpli- 
fications preserve — what we hope will be — the most essential features. 

Near the critical force, at any given time most of the interface will be 
moving very slowly if at all so that the left hand side of Eq (0) will be close 
to zero. Thus a first try might be to replace the actual dynamics with the 
adiabatic approximation that the forces at every point always balance exactly. 
Let us focus on one point r on the interface and divide the stress o~(r), Eq (||), 
into the local part Ju(r,t) and the non-local part, f a , involving w(r',t / ) for 
r' 7^ r with 

J= fdr'f dt'J{7,t'). (8) 



We then have a balance between the local force, f p + Ju and the pulling force, 

(t>{r,t) = U{r,t)+F (9) 
The adiabatic approximation corresponds to 

f p {f,t)n<t>(r,t)-Ju(r,t). (10) 
But generally, because of the non-linearities in f p (u), for a fixed <ft "applied" 
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to u(f) there can be multiple values of u which satisfy Eq (|10l); as we shall 
see these play an important role in the physics. 

At this point, it is helpful to be more concrete. Let us consider a simple 
model of the pinning consisting of pinning sites u p a {r) distributed for fixed r 
with random spacings between the v? a {r), 

T a (r)=< +1 (r)-<(r) (11) 

drawn, for each r, independently from a distribution U(T)dT. The pinning 
force f p [(r,u(r)} = except if u{f) is equal to one of the pinning positions, 
while for u(r) = u p a {f), f p can take any value between zero and a yield 
strength, f y , which is the same for each pin. A typical realization of the 
pinning force f p {u) on some segment of the interface is plotted in Fig. 2a. 
Note that for a fixed <f>, there are several possible values of u given by the 
intersection of the line 4> — Ju with f p (u). 

If cf) is increased, then the particular (history dependent) force-balanced 
position u(<f)) that the interface point is following adiabatically can become 
unstable — for example, the configuration denoted by the circle in Fig. 2a — 
and u must jump to a new position. During the jump r/^ is clearly not small 
and Eq ( |TUD will not be satisfied. But if such jumps occur more rapidly than 
the time scales of interest then their primary effect will be a time lag between 
the actual u[(f)(t)} and an adiabatic solution u a d[(fi(t)} to Eq @. A way to 
capture this feature while preserving both the physics of the time delays and 
the conceptually simplifying separation of motion into adiabatic and jump 
parts (with |j| = in the adiabatic parts for the pinning model illustrated 
in Fig H) , is to require 

u[<l>(f,t)]=u ad [(f>(f,t-t d )] (12) 

with some fixed (microscopic) delay time t d - This is illustrated in Fig ^|a. 
Note that, formally, this can be accomplished, by taking rj ^ and J(f, t) = 
J(r)5(t-t d ). 

III. Infinite-range model: mean field theory 

The above discussion in terms of the local pulling force 0(r, t) suggests that 
we could try to analyze the system crudely by assuming that the spatial 
and temporal fluctuations in <fr(r, t) are small so that <p can be replaced by 
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Figure 2: a) Simple model of the forces on one segment of an interface. 
The segment can be pinned at the positions, u p ia) of the vertical lines at 
which the pinning force can take any value up to the yield strength f y . The 
intersections of the "comb" representing the pinning force f p (uj) and the 
diagonal line ip — Jui with ip the total pulling force from the applied force 
and other segments of the interface, are the possible stationary positions of 
Ui indicated by the dots. The one of these with the smallest Ui, u™, plays a 
special role as discussed in the text. The amount Atp that <fi needs to increase 
by to depin the segment from this pinning position is WiJ. b) Dynamics of 
the same segment of the interface as the pulling force is increased. The 
actual Uiit) (dotted), the adiabatic approximation to this (solid), and the 
time delayed approximation (dashed) that is used in the analysis in the text 
are all shown. 
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some sort of time dependent average <p{t) which would then need to be de- 
termined self-consistently from the behavior of the neighboring regions that 
contribute to the stress at our chosen point r [|J . This is very analogous to 
the well known mean field approximation to conventional phase transitions: 
for example, in a magnetic system the statistical mechanics (or dynamics) 
of a single spin S(r) in the presence of a local effective field, h e ff(r) from 
its neighbors — the mean field — is analyzed and the mean field determined 
self-consistently from the condition that at all points the assumed < S > 
entering the mean field h e ff(r) = J2r> J(f — r 1 ) < S{f ) >, is the same as the 
computed <S>. One of the advantages of this approximation is that it has 
a well-defined regime of validity: in the limit that the range of the interac- 
tion is very long [or more properly that the effective number of neighbors, 
{J2f \J(f)\) 2 I Hr 1^(^)1 2 is large] then the mean field theory becomes exact. 
But — users beware! — for fixed, but finite range interactions it can still fail 
near critical points as we shall see. 

In order to obtain analytic results in at least some model, we will study 
a strictly mean field limit where each discrete segment of the interface — we 
label simply by a subscript i since it is no longer really a spatial coordinate — 
has an independent random pinning force fp(ui), of the form in Fig and 
the N segments are all coupled together by a uniform coupling 

^ = 4r ( 13 ) 



i.e., infinite range forces. (Note that Eq ( p~3|) includes a self-coupling piece 
but its effects are negligible in the desired iV — > oo limit.) Much can be done 
for general non-negative J(t) and more complicated forms of f p (u) using the 
actual dynamical evolution Eq (|3]) ||, but to keep things simple we will use 
the time-delayed adiabatic approximation discussed above with 

J(t) = J5(t-t d ); (14) 

and the form of /* of Fig ||]a with independent randomness for each i. For 
simplicity, we will focus on the strong pinning limit which corresponds to 

fy > J^rnax- (15) 

It is left to the reader to show that including some of the more "realistic" 
features within the infinite range model does not change the qualitative or 
other universal aspects of the results. 



10 



Our task is now simple, at least in principle: we assume some mean field 



<j>(t) = F + J±- y £ Ui (t-t d ), (16) 

compute the evolution of each Ui(t) (from (pi = (p for all i) from some set of 
initial conditions, and then adjust our guessed (f>(t) until the computed 

<«(*)>= 4 !>(*) 

i 

is equal to <p(t + t^) — F /J for all t. 

We first try the simplest possibility: a constant <f>(t). We can proceed 
graphically. From Fig |2|a we select for each i one of the possible stationary 
values of U{. We have many choices; the only constraint being that 

<u >= (0- F)/J. (18) 

But we must be careful: If we start choosing too many large u^s, we may find 
that < u > will become too large. We can thus ask: what are the minimum 
and maximum possible < u > for a given 0? The minimum will turn out to 
be of primary interest so we focus on this: for each i, the minimum Ui, m™(0) 
corresponds to the first pinning position — i.e. one of the {u p ia } — to the right 
of the intersection of the line f = <fi — Ju with the line f = f y that passes 
through the tips of the "comb" -representing the yield strength- in Fig. 2aQ 
Since the peaks are randomly positioned, 

< Ui >™« =< u m(ft >= (0 _ fy yj + <Wi> (19) 
with Wi the distance to the next pin which has the distribution^: 



Prob(w) 



r°° 1 


"Tn(T)c/T" 


= J w T 


T 




ensures that ' 



(20) 



i y ^ u j-max cuouico unau u, i is at a pinning position. 
The general case can be worked out similarly 

2 We use notations like "Prob(w)" to mean the probability that the continuous variable 
w is in the range w to w+dw, divided by dw; i.e. Prob(w) is the probability density (usually 
called by physicists "distribution") of w. One must remember, however, that if variables 
are changed e.g. from w to w' , then there is a Jacobian needed: Prob^') = (jrv) Prob(w). 



11 



Here the quantity in parentheses is the probability distribution that a random 
point is in an interval of width T between pins; this includes the factor of 
T/T with 



because of the presence of more points in wider intervals. Integration of 



a non-trivial result for the critical force above which no static solutions are 
possible. Note that as the interaction strength, J, is increased, the critical 
force decreases. Physically, this is a consequence of the elasticity causing the 
system to average over the randomness more effectively: pulling a stiff object 
over a rough surface is easier than pulling a flexible one. 

For F < F c the number of stable solutions, N s will be exponentially large 
with an "entropy" per segment which is of order one well below F c 
but decreases to zero at F c as most of the w^'s will then need to take their 
minimum values to ensure self-consistency. 

What happens as F is increased slowly from below F c so that the non- 
adiabaticity is negligible? Since this will certainly result in < u > and hence 
increasing, we can understand the behavior from Fig [2| As increases, 
some segments become unstable and jump to their next pinning positions. 
But this cannot happen unless they are stuck on the pin with smallest u p a , 
i.e. w™(0)- Furthermore, after a jump u will again be on the new smallest 
u p ia for the increased 0. Thus the lij's are gradually swept to their minimum 
stable positions as F is increased towards F c . 

Above F c , the segments continue to jump from one «™"(0) to the next. 
But now the time delays must play a role. If we assume a solution which 
progresses uniformly on average, < u >= vt, then 




(21) 



Eq fl20|) by parts yields < w >= T 2 /(2T) so that 




= J v t- Jvt d + F. 



(24) 
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[Note that to ensure that u does not stop between pins, we again need the 
strong pinning condition JT max < f y ] With all Ui = u™(<f)), our earlier 
analysis immediately yields self-consistency only when 



F — F r 



(25) 



v~(F- F c f (26) 



so that near criticality 
with the critical exponent 

/3 = /3mf = 1 (27) 
in this infinite range mean-field model . 



Note that a comparison of Eq ( p5|) for F > > F c and the original dynamic 
equation (13), suggests that a natural choice is td = % so that v = — for large 
F. The actual non-adiabatic processes can be seen to give rise to an effective 
td of roughly this magnitude. But the breakdown of the jump approximation 
for large F will make v not strictly linear for F > F c . Nevertheless near 
F c , the mean velocity will still be characterized by the exponent (3 = 1. A 
typical mean- field v(F) curve is shown in Fig |||. 

In addition to the critical force and the steady state velocity, in our sim- 
ple mean field model one can compute other "macroscopic" properties e.g., 
hysteresis loops as F is increased and decreased — sketched in Fig |] || — and 
responses to a time dependent additional driving force which in the pinned 
phase will depend, due to the metastability, on the past history, the sign of 
the perturbations, etc [|J. Some of the properties — like the exponent j3— 
will be universal within a broad class of mean field models, while others 
will depend on details. Nevertheless, substantial qualitative insight can be 
obtained that gives useful clues to the behavior of finite — in particular real 
two-dimensional — interfaces. Furthermore, experience with equilibrium crit- 
ical phenomena suggests that some of the asymptotic forms near F c , such as 
Eq (p5|) will be correct in realistic models if the dimension of the "interface" 



is sufficiently large or the interactions are sufficiently long range [48] [16]. We 
will return to these questions later, but for now we stick to the simple mean 
field model and ask what can be learned about "microscopic" properties, 
in particular the properties of avalanches that occur as the driving force is 
increased towards F r . 
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Figure 3: Velocity versus driving force in a typical mean field model is indi- 
cated by the solid line. Note the linear dependence of v on F just above F c . 
The dashed line is the behavior in the absence of pinning. 




Figure 4: Schematic of hysteresis loops that occur as the force is increased 
from zero to the critical force, decreased to the critical force in the opposite 
direction, and then cycled between these values. The direction of change of 
F is indicated by the arrows with the "1" denoting the first increase. 
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IV. Avalanche statistics and dynamics 

In the mean field model introduced in the previous section, the statistics and 
other properties of the avalanches of jumps that occur as the driving force 
is increased slowly can be worked out in substantial detail [59"|[H|. We will 



carry out the analysis using methods which can be generalized to provide 
useful information about the behavior with more realistic interactions. 

Let us consider what happens for F < F c when F is increased by a very 
small amount. If the increase is sufficiently small, then no segments will 
jump. But a slightly bigger increase — typically of order h — will result in 
one jump. Let us call the time of this jump t = and measure times in units 
of the delay time td so that with F held fixed after the avalanche starts, 
we have simple discrete time dynamics. The first jump can trigger n x other 
jumps at time t = 1, with these triggering n 2 further ones at t = 2 etc. 
As long as the total number of jumps is finite, then in a large system the 
mean < u >, and hence will have only advanced by an amount of order 
jj. Thus all that will matter is the distribution of segments which are very 
close to jumping, i.e. those which will jump when is increased by a small 
amount Ayv ^From Fig. 2a, we see that A<£>j = Jw^. For large N, all but 
very special ways of preparing the conditions before the avalanche starts will 
yield a distribution of these small Aipi which are independent and randomly 
distributed with (initial-condition dependent) density 

p = p(A^ = 0); (28) 

p thus measures a local susceptibiltiy to jumping. We can now immediately 
conclude something about the mean number of jumps < n t > at a time t after 
the initial jump. Since the n t -\ jumps at time t — 1 will cause an increase 
in (j) by, on average, J < T > n t -i (where we have used < T > rather than 
T since the distribution of T's for the almost unstable segments and hence 
< T > could depend on the initial conditions). This will cause, on average, 
p < T > n t -\ jumps at time t, i.e. 

< n t >= p < T > J < n t -\ > . (29) 

The crucial parameter is thus p < T > J; if this is greater than one the 
avalanche will runaway. If the system is below F c as we have assumed, it 
will eventually be stopped only when a finite fraction of the segments have 
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jumped and the system has found a stable — and more typical — configuration. 
But ifp<Y>J<l, then the expected total size of an avalanche 

s = Y, Au * ( 3 °) 

i 

is simply 

< T > 

< s >= =. (31) 

1 -p < T > J 

As we shall see however, this is not the typical size: even very close to 
criticality; i.e. 

e = l-p<T>J«l, (32) 

most avalanches will be small. 

The distribution of avalanche sizes, as well as other interesting infor- 
mation on their dynamics etc., can be obtained from generating function 
methods. Since these are a widely applicable tool, we will go through some 
of the details. For simplicity we work in units with 

< T >= 1 (33) 

and 

.1=1. (34) 

We are interested in the time evolution of the displacements, in particular 
the increments 

m t = ^[^(t) - Ui (t - 1)] (35) 

i 

— ^2 

a=l 

where the {T ta } are the magnitudes of the n t jumps that occur at time t 
(with m t — if n t — 0). The total size is simply 



£m*. (36) 



t=o 



The joint probability distribution of all the {m t } given the initial jump n = 1 
P{m t } = Prob[m , mi, m-i . . . \no — 1], (37) 
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contains the information of interest. Note that vertical bars as in Eq (|37l) 
denote "given"; i.e. conditional probability. It is useful to define a generating 
function of the distribution including all times up to T 

T 

T T {fit} =< exp(iJ2fi t m t ) > P ; (38) 
t=o 

(usually we will drop the P). Note that IV is simply the Fourier transform 
of P restricted to times < T. We can derive a recursion relation for IV in 
terms of IV- 1- For a given TOt-i, the number of jumps triggered at time T 
will be Poisson distributed with mean pmr-i i-e. 

-pm T -i 

Prob(n T |m T _i) = — (pm T _i) nT . (39) 

Tlx' 

Then, since tut depends only on tut-x, we can compute 

< e l ^ mT | {m t } t<T >=< e^ tmT \m T -i > (40) 

OO ( TIT 

= \ Pr ob(n T | m T _i) Y[ 

= exp [pm T _i (< e^ TT > -l)] 

The last equality follows from Eq. (|3"9"|); the resulting expression has similar 
mx-i dependence to the e ! '' T - imT - 1 factor in IV- 1- 
We thus find that 

r T [/i , . . . ,//t] = r T _i[/io,...,A*r-2,A T -i] (41) 

with 

A T -i = /*r-i - *p(< e iATT > -1) (42) 

where 

At = Ht- (43) 

We can now iterate with the recursion relation Eq (f42|) from the "initial" 
condition Eq fl4"3"|), eventually obtaining 

MM = T (A {M) =< e* Ao ^> T > (44) 



dT Ta Il (T Ta ) e 
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(since uq = 1). All the information has thus gone into Ao- As long as the 
system is stable, i.e. p < 1, we can simply take T — > oo to recover the full 
information. [If p > 1, then there is a non-zero (and computable) possibility 
that s = oo, and more care is needed.] 

To get the probability distribution of s, we simply set all p t — A* an d then 

Prob(s) = / e -^V A *^ (45) 

with 



dp (46) 



1 

= 2^ 

and A*(/i) the (stable) fixed point solution to Eq ( (42]) . Let us first consider 
the behavior for large s. This will be dominated by the singularity in the 
lower half complex p plane nearest to the real axis — a general property of 
Fourier transforms that can be seen by deforming the integration contour 
away from the real axis until it encounters a singularity. We only expect to 
have large avalanches for 

e=l-p (47) 

small, so the interesting regime is p and e small which suggests A* small. We 
find that to leading order, 



A* 

with 



-ie + ^-e 2 + 2ibp 



(48) 



b =< T 2 > (49) 

and the sign of the square root that has positive imaginary part for p real 
chosen. For p — > 0, this gives A* = as it must for normalization of proba- 
bility T{p t = 0} = 1. The integration contour in Eq (^) can be deformed 
so that it is dominated by the cut at p = —^ie 2 /b for large s and we thus 
have, after replacing the dummy variable p by p — ^ie 2 /b and expanding in 
small p, 

Prob(s)^ / iJlipjb e -^ e - se2/(26) (50) 



By "power counting" , we see that the branch cut must yield a -\ depen- 
dence; hence for large s, 

-se 2 /(2b) 

Prob(s) 3 (51) 

S2 
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|39|j . Note that for small e, the mean < s > is dominated by large s avalanches 
which are rare since 



Prob(s ~ — ) ~ Prob(s > -) ~ e. (52) 

We see from Eq (^l|) that these yield < s >~ - as expected. 

Two questions now arise: First, can we trust this heuristic calculation? 
And, second, is the large s, small e behavior in Eq (|51~D generic? The latter 
we have answered already: the large s behavior is the same, (except for the 
coefficient b) as long as b =< T 2 >< oo. [The reader is encouraged to find 
the behavior associated with a power law tail in the distribution of ] T[J 

As far as justifying the result Eq (0), for the case in which all the jumps 
are the same, T = 1, one can compute Prob(s) exactly by changing the 
variable of integration to z = e <A *M so that the integral in Eq ( [45]) circles an 
sth order pole at z — 0. Cauchy's theorem then yields 

Prob(s) = e- ps (P ^ S - , (53) 

with, of course, P(s = 0) = 0. From the limiting large s form 

s! sa s s e~ s V2^ (54) 

the asymptotic behavior Eq (|5T|) is found confirming the validity of the ap- 
proximations made in our first derivation of this result. Note that the p 
dependence is simply via the ^{pe~ p ) s factor in Eq (0). It is nice that the 
exact result can be found in this case, but in general, asymptotic methods 
like those we have used above give more understanding and are more widely 
applicable. Nevertheless, to convince skeptical colleagues, a few exact results 
are useful! 

In addition to the distribution of avalanche sizes, we are also interested 
in their temporal evolution. For example, one might ask what is < m t \s >, 
i.e. what is the time development of an average event of size s? This can be 
computed using the generating function. If we choose 

p T = p for r 7^ t (55) 

^Note that arbitrarily large jumps can only occur in models that also have arbitrarily 
large local yield stresses, f y . 
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and 
then 



1 dT c 



i dv t 



lit = pi + u t , 
dX 



--< n t e^ s > 



ut=0 



dv t 



(56) 
(57) 



k*=0 



whose Fourier transform in /i yields 



J m t P r ob (m t , s)dm t = J m t Prob{mt\s)dm t Vioh{s) =< m t \s > Prob(s). 

(58) 

But by multiple use of the chain rule and Eq (f42|), 



5Ao 
dv t 



' dX t -i 



x ... x 



(59) 



I'D 



evaluated with all A r = A*(/i) and all /i T = fi. 
For the constant T = 1 case we obtain, 



«9An 



du t 



v t =0 



pe 



(60) 



After shifting /i as in Eq (^) we see that, for e small and s and t large, 

1 



< m t \s >~ 



Prob(s) 



(61) 



This will be dominated by /i ~ - and hence the typical duration of an 
avalanche s will be 

t ~ y/s. (62) 



Note that this is much less then the maximum possible duration r, 
The integral in Eq (|6ll) can be done exactly (by writing fi = 



max 
x 



< mAs >~ te 2s 



8-1. 

yielding 

(63) 



for large s independent of p. Again, the behavior for large s and 1 << t << s 
is generic up to a coefficient b that should appear as in Eq. (|5T|). For the 
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particular constant T case, the exact result can be computed from Eq (|58|) 
yielding 

< m t \ s >= t + 1 - (s - forO<t<s-l. (64) 

Note that this result includes the rare tail of events which have total duration 
t of order s; the asymptotic methods needed to obtain this are a variant of 
those used above although terms neglected in Eqs (|48|J61[ ) become important 
in this region. A schematic of the evolution of a large avalanche is shown in 
Fig Note that, on average, an avalanche which is going to be large starts 
with mt growing linearly in time, independent of how big it will become; this 
is very different from the behavior of typical avalanches even near to criti- 
cality which are small. Nevertheless, typical large avalanches will have large 
fluctuations of m t away from < m t \s >; these can be studied by comput- 
ing, e.g. < m%\s > from which it can be concluded that a typical m t for a 
large s avalanche is the same order as < m t \s > except for t » The 
probability that a large avalanche stops before time t can be computed from 
Prob(m i = 0\s) which is found from f^ (fi t — > +ioo, ji T = \i for r ^ t). 



m t duration X 

max (m t )~T 
s=size ~ T 2 



scaling for L<^ 



Figure 5: Typical avalanche in a mean field model showing the increase in 
size of the avalanche, m t , at each time step. Note that the fluctuations in m t 
are larger when m t is larger. 

So far, we have not attempted to relate the local susceptibility to jumps, 
p, to the original mean field model of the previous section. In general, p will 
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depend on the past history. But on a generic approach to F c from below (e.g. 
after "training" the system by a slow increase to F c from F = —F c ), p will 
approach unity at F c and the cutoff s ~ n~p i n the jump size distribution 
will diverge. Right at F c there will be a power law distribution of avalanches 
sizes; this is analogous to the power law distribution of clusters that occur 
at the critical point for conventional percolation p0| . 

One might also hope to make connections to power law spatial structures 
at percolation |40] as well as to spatial correlations at conventional equilib- 



rium critical points. But to do this we will certainly have to move away from 
the infinite-range mean field model. This we do in the next section. 

V. Toy model: spatial structure of avalanches. 

The simple mean field model of avalanches discussed in the previous section 
can be extended in a relatively straightforward way to include spatial and 
temporal structure like that which arises for general non-negative stress trans- 
fer functions J(r,t), i.e. monotonic models. Spatial coordinates for n(f,t) 
are now needed and corresponding generating function variables p(f,t). 

In the mean field model, the probability that a segment u(r) jumps in a 
small time interval is proportional to the increase in pulling force 0(f) on it 
in that interval times a local jump susceptibility p. If we assume the same is 
true here, then we can generalize the recursion relation equation Eq (|42|) to 
include general J(f,t): 

/poo 
dr J dt'J(r -f,t' - t){exp[iT A(f , t')} - 1}. (65) 

for the case with all jump displacements T equal. Many quantities of interest 
are computable by similar techniques to those in the previous section, often 
in terms of the spatio-temporal Fourier transform of J(r,t), J(q,u). 

The mean number of jumps at r at a time t after an avalanche is triggered 
at r = at time t = is 

<m(f,t) >= f f e --V^- ±- (66) 

JqJu 1 — p I J(q, CO) 

The critical point is thus still given by 

( P TJ) mt = 1 (67) 
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with 



and the mean total size is 



j = J(g = 0,CJ = 0) 



< s >-- 



T 



1-pTJ 



(68) 



(69) 



as before. We will henceforth work in units with T = J = 1. 
A more interesting quantity is again the conditional mean: 



< n(f, t)\s > 



Prob(s) 



-ifj,s„i\* (fx) p iq-r p —iujt 



[l Jq Ju 



1 — pTJ(q, u)e 



(70) 



with A*(/z) the fixed point solution to the mean field recursion relation Eq ([42]) 
or, equivalently, Eq ( p5| ) with A and \i independent of r and t. By changing 
variables one can show that (as in the previous section) conditional statistics 
like Eq fl70|) are independent of p. 
The important physical quantity is 



K(q,w) = 1 - J(q,w) 



(71) 



which embodies the information on the space and time dependent elasticity. 
Changing variables to \i — > p + constant and noting that for large s, small p 
will dominate as before, we obtain 



< n(q, u)\s >r 



V2 



1 



-i^J2ip + .^(g, a;) 



dXe 2 ' 




A 2 



A 2 + iT 2 (g,a;) 



(72) 



For an interface with dissipative dynamics and local elasticity, in the 
absence of pinning or driving forces we have, after rescaling lengths and 
times, simply 

du o 

» = v " 



so that 



K (q, to) = —iuj + q 2 



(73) 
(74) 
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But to understand the general behavior, and to apply the results to other 
physical systems, we would like to include the possibility of long range elas- 
ticity, i.e., 

jdtJm^^ (75) 
in (^-dimensions corresponding to the static 

K s (q) ~ \q\ & = K(q,oo = 0) if a < 2 

or K s {q) ~ q 2 if a > 2. (76) 

We thus consider the general case of K s (q) ~ \q\ a with a < 2 with 

a = min(a, 2). (77) 

The total displacement a distance r from the avalanche starting point, 
during an avalanche of large size s, is obtained from Eq ( |72|) with uo = 0. We 
must thus evaluate J^e l ^' r of the last expression in Eq ([72]). For r = 0, all q 
can contribute but A is small so that we can ignore the A 2 in the denominator, 
yielding 

< Au(r = 0) |s >= T < / dt n{r = 0,t)\s >~ / . rjr . , 19 (78) 

which is of order one independent of s if 

d > d c (a) = 2a. (79) 

We thus see the appearance of a special critical dimension above which 
no segment will jump more than a few times even in an arbitrarily large 
avalanche. Indeed, we will see that above the critical dimension driven in- 
terfaces will have only bounded small scale roughness. 

For d < d c , the integral in Eq (|78|) is infinite so that small q (i.e. small 
K) dominates and more care is needed. The cutoff of / of Eq ( [72] ) when 

K ~ A yields, with a typical A ~ 4j and hence q ~ 

< Au(r = 0)\s >~ s 1 -^. (80) 

Note the appearance of a non-trivial exponent relating Am and s. It depends, 
as is usually the case for critical phenomena, on the spatial dimension. As 
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mentioned in the Introduction, Eq (|80| ) is just the kind of scaling law we 
expect near critical points. Equation (|8CiD relates characteristic scales of 
displacement to the characteristic scales of avalanche size. 

We can also say something about the spatial extent and shape of large 
avalanches, by computing < Au(r)\s >. For d > d c (a), the integral in 
Eq ( |72"D will be cutoff for q > r by the e iq ' r oscillations and hence dominated 
by q ~ - yielding 

< Au ( r )\ s >-^ (81) 
for 1 << r << where the upper cutoff arises when K(q ~ i) ~ A ~ 

For ci < d c , on the other hand, as long as r << s^, Xf a2+^| ^ or typical 
A will be dominated by g ~ and 

< Aw(r)|s >~ s 1- ^ for r << s^; (82) 

i.e. the magnitude of typical displacements is approximately independent of 
r in this range. In all dimensions, Au(r) will fall off rapidly for larger r. The 
length 

L ~ (83) 

is thus some measure of the diameter of an avalanche. 

Let us now try to interpret these results [Note that the skeptic could 
compute e.g., < [Aw(r)] 2 |s > etc. to provide further support for the picture 
below]. For d > d c , the fact that < Au(r) > is much less than unity for 
r >> 1 strongly suggests that most segments will not jump even if they are 
within r << L of the avalanche center, rather only a fraction ~ \j r d-2a, Q f 
them will jump, and these typically only once or a few times. The number 
of sites that have jumped at all within a distance R < L of the origin is of 
order R 2a << R d so that the avalanche is fractal. The total number of sites 
that jump, its "area" A is thus, by taking R ~ L, 

A~ L d f ~ s « L d (84) 

with the fractal dimension 

df = 2a for d > d c . (85) 

In lower dimensions, the picture is quite different. The approximate in- 
dependence of < Au(r)\s > of r for r << L suggests that each site in this 
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region jumps a comparable number of times ~ g 1 ^ (i /( 2a ) (with fluctuations 
around this of the same order) and hence the avalanche is not fractal but has 
area 

A ~ L d ~ (86) 
while the typical displacement is 

Au(r) ~ L c (87) 

for r < L with 

C = 2a - d. (88) 

(for d> d c , ( = 0). 

The distribution of s at the critical point that occurs at pJT = 1, is the 
same as in the mean field model. This implies that 

Prob (diameter > L) ~ (89) 

with 

K = a. (90) 
The duration of an avalanche with dissipative dynamics — corresponding to 

K(q,u) w + \q\ a (91) 

—is given simply by scaling, i.e., 

T ~ L z ~ s5 (92) 

with the dynamical critical exponent 

z = a. (93) 

Note that the relation r ~ is (not surprisingly) the same as in the infinite- 
range mean-field model. 

We have found that in our toy model, many of the properties of large 
avalanches near to the critical point (actually any large avalanche although 
they are rare away from criticality), obey scaling laws which relate various 
characteristic physical properties to each other by power law relationships. 
For example, for d < d c , an avalanche of diameter L has typical size s ~ L 2a ^ d , 
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displacement Am ~ and duration L z . This type of scaling behavior is 
one of the key aspects of critical phenomena in both equilibrium and non- 
equilibrium systems. But there is more: if we scale all lengths by a correlation 
length 

f ~ e~ 1/a (94) 

which is the diameter above which avalanches become exponentially rare, and 
correspondingly displacements by durations by £*, etc., then functions 
such as those that occur in the distribution of avalanche sizes Eq flSll), or the 



average growth of the displacements during an avalanche, < du ^^ \ s > will 
be universal functions of scaled variables such as L/£. For example, from 
Eq (|72|) we obtain, for K(q, u) = —rjiuj + Dr]\q\ a and d < d c = 2a, 

du(r,t) C^ c ( r _t s_\ 

dt |s> a 4 [? w c u t«< ) [Jb) 

with s = J dfAu(f) the total size, C u and C t non-universal (dimensionfull) 
coefficients which set the scales of the displacements and times; these depend 
on the random pinning, rj, D; etc. The universal scaling function is 



Y(R,T,m)= ! I [ dAe i4A - lUT e-^ 
Jq Jn J 



(96) 



V 2tt [A2 + (_in+|Q|«)2_ 

which depends only on the dimension, the range of interactions, and the type 
of dynamics (i.e. dissipative), as is manifested in the low frequency form of 
the stress transfer function K(q,u). As we shall see in the next section, a 
similar scaling structure is expected to exist in more realistic models. 

Let us now try applying the toy model results to the interface problem 
with d = 2 and short range elasticity, so that a = 2. This dimension is less 
than 

d short-range = ^ = 2 ) = 4, (97) 

so we have ( = 2, i.e. Au(L) » L, for large avalanches. But this is clearly 
unphysical: our original model for the interface assumed that it was close to 
flat so that, at least on large scales, we need small angles of the interface i.e. 
Vw << 1. Thus the result Eq ( p7|) violates the assumptions of our original 
model in this case. 

What has gone wrong? Is the original model bad or have we made some 
grievous errors in trying to analyze it? The answer is the latter and under- 
standing why gives some clues as to how to do better. 
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In the infinite range mean field model, the odds of any given segment 
jumping more than once are very low as long as s << \^N, since the odds of a 
specific segment jumping at all is small and each jump is almost independent 
by (justifiable) assumption. But in the finite range models, we assumed that 
this independence was still true, i.e. that the probability of a segment u(r) 
jumping in a short time interval depends only on the increase in pulling 
force, A(f)(f), on it during that interval. But this is problematic: if a segment 
has just jumped, it is much less likely to do so again until its neighbors 
have caught up. The needed increase in <j>(f) for a subsequent jump will 
typically be of order one. [Actually this is true even in our toy model, but 
the cumulative effect of many jumps causes the problem: the needed A<f>(r) 
to cause a large Au(r) that consists of many jumps should be A0(r) ~ 
JAu(r) ± 0(1) while in the toy model it is, for p = 1 where large avalanches 
can occur, A0(f) p& JAu{r) ± J^J Au{r). This difference is responsible for 
major errors when avalanches involve large A«(r)'s but not for d > d c (a) 
where Au(r)'s remain of order one.] 

Our task, then, is to somehow take into account properly the anticorre- 
lations between local susceptibilities to successive jumps. This will certainly 
involve ensuring that the statistical properties of {f p [r, u(r)]} when all u(r) 
are increased by any fixed amount are preserved; i.e. the statistical transla- 
tional invariance of the system which is lacking in our toy avalanche model. 

Remarkably, in spite of its problems the toy model correctly gives the 
statistics and properties of avalanches for large s and small e in dimensions 
greater than d c (a). The basic reason for this is the observation mentioned 
above that each segment is unlikely to jump many times during even very 
large avalanches; a real understanding, however, relies on the renormalization 
group treatment discussed briefly in the next section. 



VI. Interfaces and Scaling Laws 

Motivated by the partial success of the toy model and general scaling con- 
cepts from more conventional critical phenomena, we will now approach the 
interface problem by making a scaling Ansatz. Specifically, we conjecture that 
large avalanches near the critical point have properties which scale with their 



diameter, L, (or size s) as powers of L ||28||||2]|. In the toy model we found that 



the critical exponents which characterized these scaling laws, £, df, k and z 
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depend on the dimension of the elastic manifold, the power law decay of the 
interactions if they are long range, and on the type of dynamics, but not on 
other details of the system. This is the fundamental property of "universal- 
ity" . In contrast, the critical force and coefficients in the scaling functions 
such as Eq (j95l), will generally depend on details and hence are non-universal. 
We might thus hope that in dimensions d < d c (a) for which the toy model 
fails, the exponents will still be universal functions of d and a. 

In addition to the four exponents already introduced, there should also be 
an exponent which characterizes the correlation length: this is the diameter 

above which avalanches become unlikely (like £ ~ Smax ~ e -1 /" in the toy 
model) p7 |. As F increases towards F c with a generic past history, we 



conjecture that 

£ ~ T- 1 - . (98) 
^ (F c -F)" 1 ; 

If in the mean field model the local jump susceptibility p is smooth at F c || 
so that e ~ F c — F, then ^toy-model = ~- This turns out to be the case for 
the fuller mean field model defined in section III. 

It appears that we have five separate exponents and we would like to 
understand whether there are some relation among them. Fortunately, this 
will turn out to be the case. First, however, we can eliminate one exponent. 

If, as we expect for d < d c , ( > so that some interface segments advance 
large distances in big avalanches, then the elasticity will certainly make neigh- 
boring regions advance as well so that it seems implausible that avalanches 
would be fractal; thus we expect 

d f = d (99) 

for d < d c . 

The other relations between exponents are much more subtle. The sim- 
plest assumption — and one that is borne out — is that there is one basic length 
scale £ within the interface and a basic scale 

Au~£ c (100) 

in the direction of motion. [Note that we have been sloppy and not put 
the needed non-universal coefficients, as appeared in Eq fl95|) into Eq (|1 QQ| ) .] 
Thus, for example, if we started with a flat interface and gradually increased 
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F, when 

e = F c ~F (101) 
is small, the interface would be rough on scales less than £, with 

< [u(r) - u{f)f >~ |f - f*| 2C (102) 

for \f — f| << £, and be flat on longer scales with 

< [«(f) - u(r*)] 2 >~ e 2C (103) 

for \ f — f*| >> £ as shown in Fig. Similarly, we conjecture that for e small, 
and s large 0@] 

Prob(s) ~ -g+j-ff (s/e d/+C ) (104) 

with iY(x — > 0) — > 1 and H(x — > oo) — > — i.e., a similar form to that found 
in the toy model. Scaling implies that the exponent for the distribution of 
avalanche probability as a function of diameter in Eq fl89f) obeys 

n = B(d f + C) (105) 
We can now compute the "polarizability" of the system 

A dF] y ' 

where the arrow denotes that F is increasing. This will be given by the sum 
over avalanches which are triggered with probability pdF per segment of the 
interface yielding 

X~pf^jsH(8/?f«). (107) 
With B < 1, this is dominated by large s so that 

x „ £(l-B)(d f +C) „ e -{l-B){d !+ O r ( 108 ) 

But from the earlier results about the interface roughness, Eq ( |103| ), we also 
have 

X ~ A C C „ e -(i+<,). (109) 
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Figure 6: Schematic of a one-dimensional interface with F somewhat below 
the critical force. After starting from a flat initial configuration, the typical 
displacements of the interface are of order ^ with £ the correlation length. On 
smaller length scales, the statistics of relative displacements as a function of 
separation has a simple scaling behavior as shown. An avalanche that occurs 
as the force is increased slightly is also shown. 
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We thus have ^ 

d f + c 

and hence 

K = d f --. (Ill) 
v 

[0] Note that these relations work in the toy model for d > d c (a) but not for 
d < d c (a) due to the problems discussed earlier. 

Another relation can be derived by considering the forces one avalanche 
exerts on another section of the interface. When an avalanche that moves 
a region of diameter L ~ £ occurs, we expect it to also involve quite a few 
segments nearby to jump, but not many that are far away. In particular, 
distances R ~ £ away should be borderline. The increase in pulling force 
from the original part of the avalanche on sections ~ R away will be 

~ £ d+c /R d+a (112) 

(for short range elasticity, a = 2, the argument is more subtle). For R ~ £, 
this A0 should be comparable to the deviation e from criticality or else we 
must not have properly identified the crossover distance. Thus we should 
have 

Ap(i2~£)~^~e~r* (H3) 
yielding the scaling law 

I = <*-(. (114) 

A further relation can be derived by considering the variations in the 
"local critical forces" F^r, L) needed to make a region of diameter L advance 
by of order L^. We define 

e t (f,L)=F d (f,L)-F (115) 

to be, loosely, the deviation from "local criticality". Since the volume of the 
region through which the section of scale L will advance is V ~ L we 
expect that the random pinning forces in the region will make 



5ei(L) = ^/variance[ej(r,L)] > L(d+0/2 (116) 
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i.e. the region "does not know" F c to better than this accuracy [ I0[ . If at 
scale L ~ £, 5ei(L = £) were much larger than e, then we would expect many 
regions of size £ to be unstable so large avalanches that occurred in the past 
should have been even larger. Thus we should have 



MO < e (H7) 



yielding 



1 d + C 

- < -± (118) 

from Eq ( |116| ). This can be combined with Eq ( |114| ) to yield 

2a — d . 

C>^— • (119) 

An upper bound of 

C < 2a - d (120) 

follows from the observation that the toy model should overestimate the 
jumps of a typical segment in a large avalanche. For d > d c , we expect 
( = 0;i.e. interfaces which are flat on large scales. 

So far, we have not discussed the dynamics. One of the (helpful!) fea- 
tures of monotonic models, is that at the end of an avalanche, how much 
each segment has moved, Au(r), is independent of the dynamics although 
how long the avalanche takes, in what order jumps occur, etc. will depend 
on the dynamics @]|§. Indeed, even the exponent z can depend on the 
stress transfer J(r*, t). In particular for long range interactions whose effect 
is only felt after a time proportional to r — i.e. finite velocity of information 
propagation — one must clearly have z < 1 (this can already be seen in the 
toy model if a < 1) H. In general, however, we cannot say much about z 
without much more work. 

Nevertheless, using our scaling Ansatz, we can relate the dynamical be- 
havior in the moving phase for F just above F c to that for F < F c ||28|| P| Pi . 
In particular, we conjecture that the jerkiness of the motion occurs on length 
scales up to £ ~ ^ F J~ F y and times up to 

(121) 

while the motion is smoother on longer scales. This jerky motion will look 
like the dynamics within avalanches. If we consider, crudely, the motion to 
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be made up of avalanches of scale £ occurring at intervals of order within 
each region of diameter £, then the velocity is simply 

v „ ^ (F - F c )" (122) 

with 

= (123) 

Note again, that with the exponents of the toy model for d > d c : z = a, 
v = -, ( — and /3 = 1, the scaling law Eq ( |123| ) is obeyed. 

In the moving phase, the response of < u(r, t)—vt > to a small additional 
static force SF q e^' r , 

5 < u{q,uj = 0) > 
X® = Wq , (124) 

can be obtained exactly! [Note that here either sign of SF q is okay if it is 
sufficiently small.] If the variable change 

r ~ SF 

u(r, t) = u(f, t) + / e iq - r ^ q - (125) 

which corresponds to transforming away the response in the absence of pin- 
ning, is made, then the statistical properties of the random pinning forces 

f p [f,u(f,t)] = f p [r,u(r,t)] (126) 

as a function of u and r are identical to those of the original f p . This is a 
consequence of the underlying statistical rotational invariance of the system 
We thus find that 

< u(f, t) > SFq =< u{r, t) > 5Fq =o (127) 

so that the average response is given exactly by the second term in Eq ( |126| ), 
i.e. 

X(q) = -^—^ T7 ~ r^- (128) 

J - J(q,U) = 0) \q\ a 

Since x should scale as / L» , we again obtain the scaling law Eq ( |114j ). 
Note that since ( is defined at F c , the agreement of this result with that be- 
low threshold supports the notion that the appropriate characteristic length 
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scales above and below threshold will diverge at F c with the same exponent 
v. 

We have reduced the number of critical exponents to two basic ones; ( and 
z which relate displacement and time scales to length scales. But, so far, we 
have neither a means of computing these exponents; nor more importantly, a 
way of justifying the scaling laws beyond the hand waving arguments given 
above (or variants of these); nor even a way of understanding the claimed 
universality of the exponents. 

But with the basic scaling picture in mind, we can try and appeal to 
the framework of the renormalization group which has been so successful 
in understanding equilibrium — and some non-equilibrium — critical phenom- 
ena. There are two substantial difficulties. One is associated with the basic 
physics: we must find a way of properly dealing with two kinds of scales in 
the pinned phase near F c . First, the jumps of small segments happen on 
the basic microscopic time scale but their existence and discrete nature is 
crucial. Large avalanches last for times which scale with their size out to rg 
which is very large just below F c and thus avalanche activity spans a broad 
range of scales. But, in addition, there is the time between avalanches set by 
the rate at which F is changed; as long as this is slow enough, it does not 
really matter except that the important anti-correlations between successive 
avalanches in the same region will only be felt on this very long time scale. 

The other main difficulty is associated with the history dependence in 
the pinned phase. For example, what should one average over to get sensible 
quantities? If the stress transfer, J(f,t), is non- negative, this difficulty can 
be circumvented if F is always increasing (or always decreasing): if this is 
the case, then from any stable initial condition the pulling force on every 
segment will increase monotonically with time. This feature is important for 
the physics as well as drastically simplifying the theoretical analysis. Nev- 
ertheless, even in such monotonic models, there will be history dependence. 
But near to F c , this should, at worst, only modify non-universal coefficients 
as long as the critical force is approached monotonically starting from a much 
lower force. A natural reproducible history results from starting at F — —F c . 

For monotonic models, a perturbative renormalization group (RG) anal- 
ysis for dimensions near the upper critical dimension has been carried out 

PI |]16fl . The first result is that for d > d c (a), the decreased local sus- 
ceptibiltiy to jumping after a segment has jumped is irrelevant in the RG 
sense, except on the very long time scales during which the whole system 
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has advanced. In particular, this justifies our claim that, while the critical 
force and other non-universal properties will not be given correctly, universal 
features of avalanches such as critical exponents and scaling functions [e.g. 
H in Eq (|104|) and Y in Eq (p5|)] will be given exactly by the toy model for 
d > d c . In the moving phase, some of the effects left out of the toy model 
will be important but mean field results like (3=1 and ( = will obtain 

For d < d c many new results can be obtained from the RG analysis. In 
addition to the derivation of universality, scaling laws and perturbative com- 
putations of exponents that arise from a new critical fixed point of the RG, 
it is possible, in principle — although not yet carried out — to compute such 
quantities as the new universal scaling function that replaces Y in Eq Q9"5D, 
ant icorrelat ions between avalanches in the same region, local velocity corre- 
lations in the moving phase, etc. 

Here we just quote the results for the exponents |28| 116|. All the scaling 
laws derived heuristically above are found to be obeyed. The exponent ( is 

C ~ (129) 

to all orders in powers of d c — d; indeed, this result — which saturates the lower 
bound Eq ( |119| ) — may well be exact. Numerical computations for d = a = 1 



yield 

0.34 ±0.02 (130) 
consistent with | |8| . The dynamic exponent is found to be 

z = a _ 2ld c (a)-d\ +ou _ dn 
9 

We thus find an interesting effect: the non-linearities of the avalanche process 
cause disturbances of the interface to propagate more rapidly at long scales 
than for an unpinned interface. The velocity exponent is then, from scaling, 

a — Q 9a 

i.e. a concave downwards v(F) curve. Making the somewhat dangerous 
extrapolation to the interface with d = a = 2 , we get predictions of 

14 2 2 
z w — , C ~ -, and Bra-. (133) 
9 3 3 K J 
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At this point, except for some experiments on interfaces between two 
fluids that are being driven through porous media which give somewhat in- 
conclusive results P7|, there are not experiments, of which this author is 



aware, that test these results for interfaces. But as we shall see in the next 
section, there have been both experiments and numerical tests carried out 
on other systems. 



VII. Applications and Complications 

In the previous section we saw how scaling ideas and intuition that arise from 
the simple solvable toy model could be used to develop an understanding of 
the behavior of interfaces near to the critical driving force that makes them 
move. Renormalization group methods can then be used to carry out concrete 
calculations and justify many of the conjectures. In particular, the general 
structure and existence of scaling laws and universality follows rather directly 
from the existence of an RG critical fixed point. 

One of the advantages of this framework is that it enables us to apply 
the general structure to other systems — such as (but not limited to) dif- 
ferent dimensionalities and ranges of interactions. But in addition we can 
introduce various physical features left out of the simple models — even our 
relatively realistic Eq ([3]) — and ask whether they are relevant in the RG sense 
of changing (or destroying) the universal aspects of the critical behavior. In 
this section we will discuss several of the physical systems mentioned in the 
Introduction with an eye both to applying the ideas and seeing how they 
must be modified — or thrown away! — to account for the appropriate extra 
physics. 



A. Charge density waves 

The best and perhaps really the only real test of critical behavior near to 
a depinning transition of an elastic manifold in a random medium is that 



of charge density waves (CDW) driven by an electric field [35|. The pe- 



riodic electron density waves that occur within this class of materials are 
incommensurate with respect to the underlying crystalline periodicity in one 
direction so they could move freely in this direction — contributing to the cur- 
rent proportionally to their velocity — except for being pinned by randomly 
positioned impurities; see Fig [7|. These CDWs are thus three dimensional 
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elastic manifolds with short range interactions (the coulomb interactions are 
screened) in a three dimensional random medium. Inertial effects are negligi- 
ble. The primary difference from interfaces is in terms of the displacements 
u(r): in the frame of the crystal the random pinning forces f p [f,u(f)] are 
periodic in 

u(r) -> u(r) + Xcdw (134) 

for all r, where Xcdw is the wavelength of the CDW. Although the resulting 
behavior near the critical driving force is not much changed for d > d c = 4, 
the extra symmetry associated with this periodicity changes the universality 
class for d < d c . In particular, the exponents become @ 

C = o, 

v = —, and 
2' 

z = 2 - + 0(4 - d) 2 , (135) 
yielding scaling behavior above F c involving the velocity exponent 

f3 = zu = l-l(A-d) + 0(A-d) 2 ^l (136) 
o o 



in d = 3. Experiments |32] carried out on CDWs and numerical simulations 



yield (3 ~ 0.75 — 0.9 over 2~ decades of F — F c — surprisingly good agreement 
with the theoretical prediction. 

But should we really expect that CDWs will exhibit the critical behav- 
ior of an ideal elastic medium? In a precise sense, certainly not. Thermal 
fluctuations (or even quantum fluctuations) can cause sections of the CDW 
to overcome the barriers caused by the pinning and jump to lower energy 
local minima. For any non-zero electric field, this will cause the CDW to 
gradually creep forwards thereby contributing to the current. The critical 
behavior near the fluctuationless F c will thus be smeared out by fluctuations 
which are hence a relevant perturbation [^j, This is quite analogous to 



the role of a magnetic field near to ferromagnetic phase transitions, which 
also smears out the critical singularities. But as in the magnetic case, if 
the perturbation is small enough, critical behavior can still be observed over 
a wide range of scales and F — F c with the smearing only occurring very 
close to F c . In general, fluctuation effects appear to be quite substantial in 
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Figure 7: Schematic of a crystal with a charge density wave is incommen- 
surate with respect to the underlying lattice periodicity of the solid. The 
maxima in the difference between the actual electron charge density and that 
without the CDW are shown. An electric field pulls on the CDW which could 
move freelyexcept for being impeded by pinning forces of the impurities. 



CDW systems and the critical behavior is smeared out. But the experiments 
quoted above |32| were performed by applying an ac driving force in addition 
to the dc drive; this reduces the effects of thermal fluctuations and appear 
to yield a wide range over which they do not play much of a role. 

Another complication in CDWs — which might also be reduced by an ad- 



ditional ac drive — is defects in the CDW lattice, especially dislocations [^S 
The pinning is very weak in CDWs implying that the stresses that would 
cause dislocations to form are unlikely to occur except perhaps on long length 
scales. In this weak pinning regime, the effects of dislocations are poorly un- 
derstood. It now appears, however, that dislocations will not destroy the 



existence of the CDW phase in equilibrium [O; this is not directly relevant 



to the non-equilibrium physics of interest here but may point to progress 
also on the more relevant issues. Nevertheless, at this point, whether or not 
defects always destroy the elastic depinning critical behavior that we have 
studied is an open question. 
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B. Superconductors 

One system that is quite similar in spirit to a CDW, but for which the 
strength of the pinning forces can readily be varied , is a vortex lattice in 
a type II superconductor. Vortex lattices are pinned by impurities that im- 
pede their motion under the action of a transport current which exerts a 
force on the vortices |3| . Since the voltage is proportional to the mean vor- 
tex velocity, v(F) plots are simply volt age- current curves, F c being simply 
proportional to the critical current density. 

By making bulk superconducting samples, thin superconducting films 
or wires, or a normal layer sandwiched between two bulk superconductors, 
three, two and one dimensional systems can all be studied as well as, in the 
last two dimensional lattice of roughly parallel vortex lines with no 

dislocations allowed. In principle, many systems and regimes can thus be 
investigated, subject to the complication of non-uniform forces on the vor- 
tices, dissipative heating effects, and natural tendencies of experimentalists 
to care more about the magnitude of the critical current density than about 
what happens when the superconductor fails, i.e. when the critical current 
is exceeded! Surprisingly, although thermal fluctuation driven transitions in 
both clean and dirty superconductors have received a lot of attention re- 
cently, f34"l|12| the behavior nearer to critical currents has received far less 



If the pinning is very strong, the vortex lattice will be destroyed. What 
then (when thermal fluctuations can be neglected) will be the qualitative 
behavior near the critical force? In the case of two dimensional films in 
a perpendicular magnetic field that produces point-like vortices, the vortex 
flow just above F c will almost certainly be confined to a sparse interconnected 
network of irregular channels across the system, as sketched in Fig [8]. Some 
preliminary theoretical studies []23| and experiments of the critical behavior 
that can arise in this regime have been carried out and experiments that "see" 
the vortices performed |j25 |, but even this simple case is far from understood. 



For more complicated cases of intermediate strength pinning or when three 
dimensional effects are important, even less is known. 

One thing that is clear, both theoretically and experimentally , is 
that the critical force is history dependent. This and other history dependence- 
which can be caused by defect motion — is certain to play an important role 
in the physics. We note that in principle — and soon, if not quite now, in 
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Figure 8: A thin superconducting film in a magnetic field , H, perpendicular 
to the film. A supercurrent provides a driving force for the vortices to move 
accross the film. Just above the critical current for vortex motion, the vortices 
in most of the film can remain stationary with the vortex flow confined to a 
sparse network of vortex channels, as shown. The electric field is proportional 
to the average vortex flow rate. 
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practice |25] — measurements of local vortex flow on small scales will be pos- 
sible in some of the parameter regimes of greatest potential interest. Both 
macroscopic and microscopic information should be available in these super- 
conducting vortex systems. 

Leaving the effects of lattice defects as an intriguing puzzle, we now turn 
to another type of physical effect that we have so far left out: the role of 
inertia and other related phenomena. 

C. Cracks 

In contrast to CDWs and vortex lattices, interfaces between phases are "topo- 
logical" and hence have elasticity which is robust for weak pinning, although 
for strong pinning they can become fractal and a very different "invasion per- 
colation" behavior occurs. Various other elastic systems are also not directly 
susceptible to destruction by defects. In particular, the front of a planar 
tensile crack in a heterogeneous solid has long range \ elastic forces medi- 
ated by the elasticity of the solid which act to keep the crackfront roughly 
straight; see Fig|] || [50], But random variations in the local toughness— 



the energy per unit area needed for crack growth — act to deform the crack 
front. This system is thus an example of an elastic manifold with d = 1 and 
a = 1. 

Numerical studies PH [|3*1"|| with quasistatic (i.e. instantaneous on the time 
scales of the crack motion) stress transfer which corresponds to J(r, t)a5(t) / r 2 
and hence is a monotonic system, yields results for (, z, ,u, and f3 in excel- 
lent agreement with the d = 2 — e expansion results |16 |. But the resulting 



z ~ | is less than one || [HJ. Therefore, even if the basic propagation of 
disturbances along the crack front is slow on small scales (so that the qua- 
sistatic approximations would seem to be justified), near the critical point 
disturbances would propagate arbitrarily fast since the characteristic scale 
of the "velocity" along the crack front would scale as which diverges as 
£ — > oo. This is clearly unphysical and so elastic wave propagation effects 
must change the asymptotic critical behavior. The simplest way to take the 
time delays associated with wave propagation into account is to use 

J(r,t)~\s(t-r/c) (137) 

with c an elastic wave velocity. This is still a monotonic model, so only 
the dynamic exponent z will change. Theoretical arguments and numerical 
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Figure 9: A solid with a crack under tensile loading. If the load is increased, 
the crack front can progress through the solid. In some circumstances, the 
crack is confined to a plane; this is the case discussed in the text. 
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simulations M for this case yield z — 1 exactly — so there is no problem with 



causality. From Eq (|123|) , we obtain (3 = 1. 

But real elastic wave propagation is much more complicated If 
the crack grows by a segment of the crack front jumping forward only to 
be stopped by a tougher region, then a point on the crackfront a distance 
r away will initially feel no change in stress. Furthemore, when the lon- 
gitudinal waves with velocity q arrive, the stress that tends to open the 
crack — analogous to <j(r,t) for the interface model — will actually decrease 
since J(r, t > t/q) < [0, HQH]! Only when the Rayleigh waves with 
slower velocity cr < c t < q arrive will the opening stress become positive. 
But when this occurs, there will be a rapid jump to a large peak value of 
the opening stress followed by a gradual fall off to the eventual static stress 
increase which is 

/•DC | 

J( r ,u = 0)= J(r,t)dt~—. (138) 
Jo H 

This is highly non-monotonic behavior! Such stress overshoots can cause seg- 
ments of the crack front to jump forward that would not have been triggered 
by the static stress changes. How often this occurs will depend on the heights 
of the stress peaks relative to the static stress increases — i.e. to the size of 
the overshoots. If the jumps of segments of the crack are slow and smooth 
(with duration long compared to the sound travel time across the jumping 
segment) then the stress overshoots will only occur far from the jumping 
segment and have small amplitude. Their effects should then be small and 
the quasistatic behavior should be observable — except, as we shall see, very 
near F c . But any stress overshoots will at least occasionally cause some extra 
jumping and we must understand their effect. 

Some intuition can be gleaned by considering what would happen if the 
stress increases never decayed to their static value, but stayed at their peak 
strength. This would correspond, roughly, to increasing the elastic inter- 
actions of the crack with itself. As we saw in the mean field model, such 
enhanced elasticity causes the critical force to decrease because of more ef- 
fective averaging over the randomness [see Eq fl20D]. Thus we would expect 
avalanches to run away at a force below the quasistatic F c . A careful analysis 
of the effects of realistic stress transfer peaks followed by decay towards the 
static stress shows that the same basic picture obtains J5|]||. 

Consider a segment which, at the end of the quasistatic approximation 
to an avalanche, has a final stress on it which is very close to being enough 
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to make it jump again. Then any overshoot in the stress on this segment 
that occured during the actual avalanche will, if the resulting peak stress is 
above the quasistatic final stress, cause it to undergo an extra jump. If, as 
is indeed the case, there is a small density of these extra jumps roughly ran- 
domly distributed within the avalanche, each of the extra jumps will trigger 
extra avalanches, for which one again has to consider the effects of the stress 
overshoots which make more avalanches, etc. One can show that this process 
will always run away for a sufficiently large avalanche — independent of how 
small the stress overshoots might be || H . 

The stress peaks are thus a relevant perturbation which will change the 
critical behavior. Indeed, the runaway of large avalanches implies the coexis- 
tence of moving and stationary states and strongly suggests that the depin- 
ning transition in the presence of stress pulses will be either first order with 
a discontinuous and probably hysteretic v(F) or have a continuous v(F) but 
with a smaller critical force and different critical behavior ||. 

What happens if the stress overshoots are small, say with strength pa- 
rameterized by S? Then only avalanches with 

L > L S ry- (139) 

will typically be much affected by the stress peaks. The exponent y% is the 
RG eigenvalue for the relevant perturbation, £; it has been computed analyt- 
ically || H] and checked numerically |§j for some types of stress overshoots, 
but it's value is not yet understood analytically for the type of long-tail stress 
peaks that occur in the crack problem H, [U 



Experiments on very slowly advancing cracks confined to a plane |17| yield 
an estimate of the roughness exponent which characterizes the deviations of 
the crack front from a straight line, of 

(crack front « 0.55 ± .05 (140) 

substantially larger then the quasistatic prediction of |. However the large 
corrections to scaling predicted by the RG analysis HHjlJ may account for 
this discrepancy; beware of quoted error bars for exponents! It is possible, 
however, that real elastodynamic effects play a role even for slowly moving 
planar cracks. 

A related challenge is to understand the crack front distortions which 
manifest themselves in the roughness of the fracture surface left behind af- 
ter a crack front passes: i.e. after the material is broken |18|l HI • Again, 
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elastodynamic effects may play a crucial role. In certain situations, cracks 
can form multiple branches -such as shattering of glass-, or just a few small 



branches near the crack front \T5\. Under what circumstances crack branch- 
ing effect might affect the onset of crack motion or the large scale shape of 
fracture surfaces is another open question. 

D. Faults 

So far — except for crack front roughness — we have only discussed systems for 
which the simple measurements are of macroscopic properties such as v(F). 
But there is one natural system in which, even though the length scales are 
huge, the measurements are of "microscopic" type: specifically, the statistics 
and other properties of "avalanches" . 

A crude model of a geological fault is the motion of two large blocks of 
crust in contact along a disordered but roughly planar surface — the fault— 
which move relative to each other; see Fig [HJ The motion is driven by forces 
exerted from far away (e.g. from the viscoelastic region beneath the crust) 
that are transmitted to the fault by the elasticity of the blocks. Rather 
than being driven by a constant force, the blocks are driven at a fixed time- 
averaged velocity that is extremely slow — of order millimeters to centime- 
ters per year — compared to other characteristic velocities — e.g. the speed of 
sound in the rock ^|. Segments of the two dimensional fault plane interact 
with each other by long range 4j stress transfer -as can be seen by dimen- 
sional analysis- and are pinned by heterogeneities on the two fault surfaces 
that have to rub past each other Q. 

Earthquakes, of course, are just the "avalanches" that occur when some 
segment of the fault becomes unpinned and jumps forward triggering others 
47], [U| . The range of length scales of slipping regions is huge: from meters to 



hundreds of kilometers, with slips — changes in relative displacement across 
the fault plane, i.e. Am — from millimeters to tens of meters. 

Most theoretical approaches to modelling faults have involved physics 
driven by friction laws and inertial effects [ 20[| , with intrinsic heterogeneities 



playing a secondary role if included at all. An alternative approach [||], 
motivated by the systems we have discussed here, might be to start from the 
opposite end, a strongly disordered fault with quasistatic stress transfer, and 
then bring in other features such as elastic wave propagation (which is the 
manifestation of inertia) and frictional weakening. This has the advantage of 
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Figure 10: Schematic of two segments of crust which move relative to each 
other along a fault plane in a sequence of earthquakes. The driving forces are 
transmitted to the fault by the two halves of the crust moving with relative 
velocity, v. 

building on an established theoretical framework in which the importance — 
or lack thereof — of various features can be considered. This is particularly 
useful as the parameter space in even relatively simple models is very large 
and thus hard to explore numerically — at least in ways that might convince 
a skeptic of the interpretations or predictions! 

If the forces that one side of the fault exerts on the other are independent 
of the history and the local slip velocities dU g t '^ and the stress transfer is 
quasistatic, i.e. J qs ~ 5(t)/r 3 , then this system falls into the class of gener- 
alized interface models with d = 2 and a — 1. But for this a, two is exactly 
the upper critical dimension so we can use the results of the toy model. The 
fault system acts as if it were driven just below F c (by an amount that goes 
to zero for large system size) so that we can set e = |i9|]. [Note that there 
should be logarithmic corrections to various quantities, as is usual at critical 
dimensions, but their effects are small]. The toy model is sufficiently simple, 
that one can include a realistic way of driving the faults, but for now we 
consider the simpler idealized infinite system. 

Several results can immediately be used ||: the area of the region of 
the fault that slips in an earthquake will scale as the square of its diameter, 
L, the typical slip in this region will depend only logarithmically on L (i.e. 
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( = 0), and the size, in this context called the moment, will scale as 



M = J drAu(r) ~ L 2 . (141) 

Note that the moment magnitude often quoted in newspapers is 

m = log 30 M + constant, (142) 

the base 30 being for historical reasons, particularly consistency with the 
older Richter magnitude scale. The probability of different sized events in 
the simple model is 

Prob (moment > M) ~ — — (143) 

with B = \. 

How do these predictions compare with observations? The moment of 
quakes is the best measured quantity, although for medium and larger size 
events, the duration r can also be measured well. The diameter — or generally 
the dimensions of the slipped region — and the slip Am can only be measured 
directly if the quake involves slip where the fault breaks through the surface 
of the earth. However, earthquake data are usually interpreted in terms of a 
crack picture in which £ = 1 and z = 1 so that 

M crack ~ L 3 ~ r 3 (144) 



J3"]. This scaling appears to be reasonably well justified observationally for 
large earthquakes and the M ~ r 3 scaling perhaps also for intermediate size 
events. But for small events — i.e. most of them — only the moment can be 
measured reasonably reliably. Thus, perhaps, our M ~ L 2 ~ r 2 might not 
be ruled out for small events although it probably is for large ones. Note 
however, that in quakes of magnitude m = 7 or so and above, the linear 
size (e.g. the diameter) is comparable to the depth of the crust (and other 
relevant length scales) so different scaling laws should in any case be expected 
for large earthquakes with a crossover between the two regimes [f|2| . 

The question of earthquake magnitude statistics is a complicated and 
controversial one. The famous Gutenberg-Richter law states that the 
distribution of all earthquakes approximately satisfies Eq ( |143| ) with a i? ~ |, 
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although it changes somewhat — perhaps for the reasons mentioned above — 
around magnitude seven. The data cover over twelve orders of magnitude in 
the moment, equivalent to, assuming the crack scaling Eq ( |144 ), four orders 
of magnitude in length scale. Understanding the Gutenberg-Richter law is 
a very interesting problem, tied closely to that of understanding the appar- 
ent fractal distribution of some geological features such as fault networks 
|n| . But our subject here is individual faults — at least if they really exist 
as well defined entities, a question which is also controversial. Observations 
of some highly disordered faults have been fitted — over a reasonable range — 



with B ~ 0.5 — 0.6 ^TJ. This is perhaps encouraging although the apparent 
rough agreement with our Eq (|143|) may well be fortuituous. Most faults, 
however, exhibit quite different behavior: a small regime (if at all) of power 
law statistics that can be fit by somewhat larger B's, then a gap with few 
events and a narrow peak (which dominates the total slip) at a characteristic 
earthquake size in which the whole fault section slips ||21||. This behavior, 



sketched in Fig |TT]b, appears to be very different from the pwer-law scal- 
ing behavior we have found in our simple quasistatic heterogeneous model. 
However, behavior qualitatively similar to Fig [TT| b — with large characteris- 
tic earthquakes — has been found in simulations of models with inertia and 
frictional weakening but no intrinsic randomness although analytic under- 



standing of these results is rather limited p0 



Can we understand qualitatively how both power-law and characteristic- 
earthquake types of behavior might arise starting from our simple randomness- 
dominated picture? Elastodynamic effects, as for the crack front problem 
discussed earlier, will result in peaks in the dynamic stress transfer that are 
larger than the static stress transfer. In addition, frictional weakening — 
the tendency for dynamic frictional forces to be less than static frictional 
forces — will also be present. Both of these effects will cause extra segments 
to slip that would not have slipped in quasistatic events. These, as for the 
crack front, will cause runaway of large earthquakes which will eventually 
be stopped only by strong "pinning" caused by the boundaries of the fault 
section, or by the unloading of the shear stress that is driving the fault. If 
the stress overshoot and frictional weakening effects are small, there should 
be a wide regime of power law scaling with B = | which can extend out to 
the largest quakes in the fault section. But if these effects are strong — as 
one might guess would be the case for a weakly disordered fault that is close 
to planar — intermediate size events will not occur and the behavior will be 
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Figure 11: Schematic of two types of earthquake statistics observed on 
different geological faults; see [20] . The probability of quakes with a given 
moment is plotted on a log-log scale, a) A fault with power law statistics of 
events, b) A fault that exhibits "characteristic earthquake" behavior which 
refers to the peak in the distribution for large events of a characteristic size, 
with not many intermediate size events occuring. 
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qualitatively similar to the characteristic- quake behavior observed in many 
faults ||. Which behavior obtains would be determined by how the length 
scale above which events typically runaway compares with the length of the 
fault section. 

Although this basic scenario of fault dynamics that was found in the sim- 
ple heterogeneous fault model may have nothing to do with what happens in 
the real earth, analog laboratory systems might be found — or synthesized — in 
which some of these ideas could be tested. 

In these lecture notes, we have outined a framework for studying the 
non-equilibrium " critical" behavior that occurs near the onset of macroscopic 
motion in many driven systems. This enables us to understand the origins, 
nature, and statistics of avalanche-like events that can occur in these systems, 
as well as other qualitative and quantitative aspects of the critcal behavior. 
But the examples given in this section have also illustrated a key point of 
these lectures: to be really useful, a phenomenological framework should be 
broad enough and robust enough to show the roots of its own failure. With 
judicious choice of which experimental systems to focus on - and some luck 
- this should enable enough predictions to be made that theoretical results 
and underlying assumptions built into models are falsifiable! 
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